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Abstract 


Pd/n-Si Schottky contacts have been fabricated by depositing a 
palladium film of ~ 80 nm thickness in vacuum ~ 10~ 6 m r , using an elec- 
tron beam evaporation source onto a cleaned n-type <111> silicon wafer 
and subsequent vacuum annealing at 200°C~500°C for different lengths of 
time. Their current- voltage (I-V) characteristics have been measured in the 
tomperatuic range 80K-300K and analysed on the basis of the thermionic 
emission-diffusion(TED) mechanism and parameters, viz., barrier height , 
ideality factor , saturation current etc. are determined. The zero bias 
barrier height (<^ 0 ) is found to decrease while ideality factor^) increases 
with decrease in temperature. Further , the changes in ^ and rj are quite 
significant at low temperatures. This abnormal behaviour has been ex- 
plained by invoking the concept of barrier inhomogeneities and describing 
them by double Gaussian distribution of barrier heights with parameters as 
: zero bias mean barrier height (<^ D ) 0.93V and 0.63V, standard deviation 
(a 0 ) 11.2 and 67 6mV in the temperature range 165-300K and 80-165K, 
respectively. 

It is shown that thermal annealing at 400°C , 450°C and 500°C of 
Pd/n-Si invariably leads to the formation of silicide of composition Pd 2 Si 
possessing hexagonal structure (Fe 2 P type ) with parameters a=13.055A , 
c=27.49A Also the growth occurs preferentially in certain directions such 
that (11.4) or (00.1) planes of hexagonal Pd 2 Si lie parallel to the underly- 
ing (111) Si substrate in the beginning, but, shows polycrystallinity with 
time at all the temperatures. The reflectance is shown to usually decrease 
with the growth of Pd 2 Si. In the intermediate stages a slight increase in 
reflectance observed in particular cases is attributed to the disposition of 
Pd 2 Si preferentially on the underlying (111) Si substrate , As the time 
lapses , Pd 2 Si layer exhibits polycrystallinity at all the temperatures and 



corresponding reflectance decreases - taking a constant value eventually. 

Finally, simulation studies have demonstrated the dominance of 
TED mechanism for current transport through the Pd/n-Si Schottky con- 
tacts. Also, negative temperature coefficients of mean (<£(*,) and variance 
(crl) of the Gaussian distribution function (representing the barriei inho- 
mogeneities at the Schottky contact) are shown to have opposite effects on 
the temperature dependence of apparent barrier height . While the former 
causes decrease in barrier height , the later compensates the change by 
a fixed amount. The overall effect is that Gaussian distribution function 
shifts towards lower barrier height and at the same time become narrower 
with rise in temperature, i.e., effective reduction in barrier inhomogeneities 
occurs. 
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Chapter 1 

The Schottky Barrier 

1.1 Introduction 

Within the framework of solid state and interface physics , the solid -solid inter- 
face is traditionally studied the most . One reason is the enormous importance of 
metal-semiconductor (MS) junctions and semiconductor hetero junctions for device 
physics. The metal-semiconductor junction , also known as Schottky diode , is an 
active device capable of giving a strong non-linear response under appropriate con- 
ditions [1-2]. It is found in many metal semiconductor contacts and the underlying 
physics forms the basis for their rectifying properties.There have been three major 
periods of industrial importance for the metal-semiconductor contacts - in early 
twentith century they were used as radio detectors , then as radar detectors during 
the second world war and later on as clamping diodes, microwave diodes and gates 
of microwave transistors [3]. 

The first metal-semiconductor device was made, during the early days of teleg- 
raphy , by pressing a sharp metallic wire to an exposed semiconductor surface, for 
use as radio- wave detector [3, -S J. But, they proved to be highly unreliable in 
their characteristics and were subsequently replaced by reproducible planar metal- 
semiconductor interface , produced by depositing films under high vacuum , with 
an uniform contact potential and current distribution over the junction area. 

1.2 The formation of a Schottky barrier 

An ideal surface essentially is one having atomically sharp profile with no in- 
termixing of the components on the either side , e.g., metal and semiconductor in 


1 



The Schottky Barrier 


2 


case of a MS contact. Its electronic properties are described primarily by the work 
function (j) m of metal , electron affinity \s of semiconductor, position of their Fermi 
levels and band structures. When the two are brought into contact , equilibrium re- 
quires that the chemical potential on both the sides of interface should be the same. 

This means that the Fermi levels in the metal and the semiconductor be aligned. The 
matching of the Fermi levels requires charge to flow from one side to the other In 
the process, free electron concentration in the semiconductor near the boundary 
region decreases and, in turn, causes corresponding shift of the Fermi level. Since, 
the seperation between the conduction band edge and Fermi level increases 
with decrease in electron concentration and Er becomes the same throughout at 
equilibrium, the conduction band edge Ec bends upward (Fig. l.l(i)). If the band 
gap of the semiconductor is not changed, the valence band edge will also move up 
parallel to the conduction band edge Ec- Also, the electron affinity (x s ) of the 
semiconductor is assumed to remain unchanged. Hence, the vacuum level should 
similarly bend upwards to maintain the value of Xs- Thus, the vacuum level remains 
continuous across the transition region. Depending upon the work function of the 
metal and electron affinity of the semiconductor different situations may arise as 
shown in Fig. 1.1 (ii) 

A dipole layer is built up at the interface due to the flow of the charge ,the 
size of which can be estimated by considering the screened coulomb potential of a 
point charge at a distance ’r’[6,7], 

<f>(r) = (g/r)exp(—r/r TF ) (1.1) 

where q is the electronic charge and ttf = 0.5(n/a o ) _1 ' /6 is the Thomas Fermi 
screening length; a 0 is the Bohr radius and n is the electron density. With n~ 
10 22 cm~ 3 (i.e., for metals), t F f is ~ 0.55 A and with n~ 10 18 cm~ 3 (e.g. ,for semiconductor) 
is ~ 5 00 A . Thus a normal space charge layer (depletion or accumulation type) is 
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Figure 1.1: Schematic diagram of band bending before and after metal- semicon- 
ductor contact (a) High work function metal and n-type semiconductor ,i-e., <j> m 
> <t> s ,( b) low work function metal and n-type semiconductor , (c) high work func- 
tion metal and p-type semiconductor , (d) low work function metal and p-type 
semiconductor [ 6 ] 
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formed. The amount, of band bending in the depletion region is correlated with the 
positive space charge of ionized donors via Poisson’s equation 

drV _ -p(z) 
dz 2 eeo 

whexc V is the potential, p(z) is the space charge density at a distance z from the 
interface. The positive charge is balanced on the metal side by a corresponding 
excess electronic -charge extending over a subatomic distance. The maximum band 
bending eV& is just equal to thp difference between two vacuum levels(see, e.g. Fig 
l.l(i)).Thus, 


vVfc = 4>m- <f>s (1.2) 

where is the potential barrier encountered by an electron moving from the semi- 
conductor to the metal. 

The potential barrier efo which has to be overcome when an electron is excited 
from thp metal into the conduction band of the semiconductor is given by 

4>b = <j>m ~ ^ (1-3) 

Using equation 1.1 and substituting X, = 4> s — <f> n . we have 

r > 


4>b = qVb + <f>n (1.4) 

Schottky assumed the semiconductor to be uniformly doped , giving a 
constant charge density in the depletion region. The electric field strength therefore 
rises linearly with distance from the edge of the space charge layer [4], Mott assumed 
a thin layer devoid of charges sandwiched between a uniformly doped semiconductor 
and the metal. The electric field strength within the so called depletion region is 
constant [4]. The depletion region has few mobile carriers and hence has very high 
resistance . Consequently , a large part of the externally applied voltage drops across 
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the depletion region. If a negative voltage (V f) is applied at the semiconductor , the 
depletion region width decreases and voltage across it falls from Vf, to Vfe-V/.The 
electrons now see a reduced barrier and move from the semiconductor towards the 
metal causing a current flow . Since no voltage drop occurs w’thin the metal , 
<j>b remains unaffected and there is no electron flow from metal to semiconductor . 
Thus MS junction exhibits rectifying properties. 

The difficulties involved with this simple approach are obvious. Firstly, 
the barrier height in practise does not increase linearly with <j) m in accordance with 
> equation (l.'\) and (l.#*). Secondly, the real work function is not a constant but 
depends on the surface orientation and gets strongly affected by types of rearrange- 
ment of surface atoms , e.g. , relaxations and reconstructions^] . Here the nature 
of surface and interface comes into picture. 

1.3 Surface and Interface states 

The bulk crystal structure is broken at the surface so that the coordination 
of species becomes markedly different. For the topmost surface atoms the bonding 
partners on one side are missing in total .Therefore , their wave functions have less 
overlap with that of the neighbouring atoms.The splitting and shift of the atomic 
energy levels is thus smaller at the surface than in the bulk. According to Bardeen 
broken (or dangling) bonds at the semiconductor surface give rise to localised energy 
states called surface states. These states are distributed in the band gap and are 
characterized by a neutral level <j> 0 . That means when there is no band bending in 
the semiconductor, the surface states are occupied right upto <b 0 making the surface 
neutral , i.e., a flat band condition (Fig. 1.2 )[^] . The states below <p 0 are donor 
like while those above behave like acceptors . Depending upon the type of surface 
states and on the position of the Fermi level , the surface state may carry a charge 
, which is screened by an opposite charge inside the semiconductor .If an example of 
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Figure 1.2: (a) Qualitative explanation of the origin of surface states in the tight, 
binding picture . Depending on their origin , these states have acceptor or donor 
like charging character (b) Electron energy band diagram of n-type semiconductor 
with surface states (i) under flat band ,(ii) in thermal equilibrium with the bulk, 
and (c) in contact with metal [4,&J 
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a n-type semiconductor is taken , we have acceptor surface states ; their energetic 
positions being fixed with respect to E c . Flat bands upto the surface would position 
the Fermi level far above the acceptor type surface states; filling them completely 
by electrons and building up an uncompensated charge density.This is energetically 
unfavourable and hence band bending is there at the clean surface even in the 
absence of metal contact [6 j . 

A different approach is via the mathematical treatment of Schokley [ g ] . Here, 
a monoatomic linear chain terminated at one end is considered and the Schroedinger 
wave equation is solved for an electron under nearly free electron approximation. 
Possible solutions , with real electronic wave functions , are standing Bloch waves 
within the chain , matched to the exponentially decaying tails on the terminated 
side.Hence, the electronic energy levels are slightly modified from those of infinite 
bulk. Additional surface solutions become possible if complex wave vectors are 
allowed. The calculation then yields one single electronic surface state located in 
the gap of bulk states.Electrons in these states are localized within a couple of 
Angstroms of the surface plane.The results can be generalised to the 2D surface[6] 

Numerous experiments on the wide variety of MS systems prepared under 
UHV conditions indicate that the deposition of metal films produce interface states 
which determine the position of the Fermi level at the contact [6] . According to 
Heine [=C] metal (Bloch) wavefunctions tail into the semiconductor in the energy 
range in which the conduction band of metal overlaps the forbidden gap of the 
semiconductor. Fig.l.3(b) represents the situation in greater detail. The break- 
down of the perio deity at the interface introduces exponentially decaying interface 
states with imaginary wave vectors, which fills the energy gap of the semiconductor 
s ymm etrically with respect to the band edges E c and E t , . In a simple ID model 
Hieir dis^rsion curve E -({ (Fig. 1.3(a)) liidicates the range of theoretically possible 
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Figure 1.3: Origin of metal induced gap states at MS interface (a) Dispersion curve 
of a ID model of semiconductor ; the real band structure with E c (k) and E„(k) as 
conduction and valence bands respectively. D t , s is density of the virtual states (b) 
Tailing of metal Bloch state into semiconductor 

states known as "Virtual induced gap states" (VIGS). In the case of free surface , 
only a single energy level exists , but in the present case ,the Bloch states of metal 
conduction band have to match resulting in a broad continuum of states in the en- 
ergy gap of the semiconductor. Or else, the continuum of metal states leaks into the 
VIGS which are themselves derived from the semiconductor band structure . Which 
of these states actually exist depends on the boundary conditions at the interface. 

Fig.1.4 shows barrier heights as predicted by the Schottky , Bardeen models 
together with experimental findings for various metal-silicon contacts. Accordingly , 
the experimental result lie in between the two extreme approaches namely Schottky 
with no interface states and the Bardeen with high density of interface states. If 
the density of surface states is high , the charge exchange takes largely between 
the metal and the surface states and the space in the semiconductor nearly remains 
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Figure 1.4: Barrier heights of Si-Schottky contacts versus metal work function. For 
comparison the predictions of the Schottky (no interface states) and the Bardeen 
(high density of interface states) model are given [6j. 

unchanged. As a result , barrier height becomes independent of the metal work 
function and is given by [4 ] 


<t> b = E g - <j> 0 (1.5) 

Here the barrier height is said to be pinned by surface states and equation(1.4) 
is known as the Bardeen limit. The assumptions of Bardeen are,of course, incorrect , 
since the surface states of the clean surface are strongly affected by metal deposition. 
Taking into account the interface states , the Schottky barrier diode can be described 
as follows (Fig. 1.5): 

When the metal comes into close contact with the semiconductor surface , 
they’ll form chemical bonds . The distribution of intrinsic surface states of the clean 
semiconductor surface will be changed . Charge will flow from one side to other 
due to formation of bonds which is described by the formation of a dipole layer 
of atomic dimensions. The direction of the dipoles is determined by the difference 
between the work function and electron affinity. Formation of new interface states 
is expected in all the cases. In a more extensive band scheme of the MS contact , 
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Figure 1.5: Band diagram of a metal-semiconductor junction in which the interface 
region (width ~5 A) is taken into account explicitly 

interface region could be included with its atomic dimensions. Using the metal work 
function(^ m ) and the electron affinity (% s )of the semiconductor , dipole energy (6) 
can be associated to the interface layer . The total charge balance at the interface 
now includes the space charge in the semiconductor and the charges present in the 
metal and those located in the interface states. In thermal equilibrium , no net 
current flows . It is obvious that the band bending and in turn , the Schottky barrier 
height depends on the nature of the interface states.The origin and character of 
these states need to be further investigated for better understanding of the Schottky 
barrier formation. 

1.4 Conduction Mechanisms 

The various ways in which electrons can be transported across the MS junction 
under forward bias is shown schematically in Fig.1.6 for an n-type semiconductor. 
The inverse processes occur under the reverse bias condition . The mechanisms in- 
volve: (a) emission of electrons from the semiconductor over the top of the barrier 
into the metal (b) quantum mechanical tunneling through the barrier (c) recom- 
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Figure 1.6: Transport processes in a forward-biased Schottky contact [4]. 

bination in the space charge region (d) recombination in the neutral region(hole 
injection) [4] 

1.4.1 Thermionic emission-diffusion Mechanism 

The motion of an electron crossing over the barrier from the semiconductor 
into the metal , moving through the high field depletion region, is governed by the 
drift and diffusion processes [ U']. Density of the available states, in fact controls the 
jumping (or emission ) of electron into the metal. The two processes , the emission 
over the barrier and the drift and diffusion in the depletion region are effectively 
in series and the one which offers the higher resistance determines the current. In 
Bethe’s TE theory [7j , it is assumed that the barrier height is much larger than 

kT and the effect of drift and diffusion in the depletion region is negligible. Hence 
, the current flow depends solely on barrier height. Schottky in his diffusion 
theory had also included the effect of electron collisions within the depletion region 
and assumed the carrier concentration at any distance to be unaffected by current 
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flow.Both the approaches led to the similar current expressions. According to the 
synthesis of the two theories proposed by Crowell and Sze [. > ] the net current due to 
flow of electrons from semiconductor to metal by overcoming the potential barrier 
of height (fib at a forward bias V is given by [4,7] 

J = A i A"T 3 exp(^)[exp(^)-l] (1.6) 

where A** is the effective Richardson constant, Ad is the diode area , T is the tem- 
perature in Kelvin , k is the Boltzmann constant and q is the electronic charge.For 
free electrons Richardson constant is 1.12 x 10 6 Am~ 2 K~ 2 . However , the barrier 
height is known to increase with forward bias. If the variation is assumed to be 
linear , we can write [7] 


MV) = <t>bo + lV (1.7) 

where (fibo is the barrier height at zero bias and 7 (= is positive. Substituting 
the value of cfib(V) into eqn (1.6) , we get 

I = I s exp(^)[l - exp(^^-)] ( 1 . 8 ) 

where 

I, = A d A"T 2 exp(^^) (1.9) 

and 

i = i-. 

The parameter 77 is called the ideality factor and usually has a value greater 

than unity. Further , the presence of a series resistance(R s ) causes significant voltage 

?- 

drop at large forward currents.As a consequence , the current equation (l.£) gets 
modified to 
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Figure 1.7: Field and thermionic-field emission under forward bias 


r x ,g(V — /i2 s ). r , — ^(F — 7i? s ) 

I = I s exp( f]kT — ex P{ )] * (1.10) 

Also, the ln(I) vs V plot begins to deviate from linearity and exhibits saturation. 


1.4.2 Tunneling through the Barrier 

This mechanism prevails when the semiconductor is heavily doped such that the 
Fermi level lies above the bottom of the conduction band and depletion width is 
very small. It may be possible for electrons with energies below the top of barrier to 
penetrate by the well known quantum-mechanical tunneling. In the case of degen- 
erate (heavily doped) semiconductor at low temperatures , the current in forward 
direction arises from the tunneling of electrons with energies close to Fermi energy 
in the semiconductor . This process is known as field emission (FE) . At raised 
temperatures , electrons excited to higher energies have higher tunneling probabil- 
ity as they experience thinner and lower barrier . This is known as thermionic field 
emission (TFE) [4] (Fig. 1.7). If the temperature is still raised further , a point is 
reached in which virtually all the electrons have enough energy to go over the top 
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of barrier; it is pure thermionic emission with negligible tunneling.The tunneling 
current through the barrier is of the form [4, 1 , ?'t] 



E a = E o0 coth ( |^) (1.12) 

and l t0 is the saturation tunneling current. E ao is defined as 

= A =18.5X10“ 15 t /— eV (1.13) 

47 r V rii e *e s V m r € r 

where m e *(=m r m 0 )is the effective mass of the electron and e s (=e r e 0 ) is the 
permittivity of the semiconductor, m 0 is the electron rest mass and N d is the donor 
concentration in m~ 3 . The process is FE when E && >>kT and is TFE at E « kT 
. TFE also causes a reduction in barrier height by [4] 

A <p = (?) 3 eJvJ (1.14) 

II 

where stands for the voltage corresponding to the band bending. Eqn (1.1 1 j 
suggests that In I t vs V plot should yield a straight line with slope giving E 0 and 
intercept at zero bias the saturation current l io . 

1.4.3 Carrier Generation and Recombination 

At zero bias there exists a thermal equillibrium , as the rate of electron-hole pair 
generation in the depletion region is balanced by the rate of recombination. With 
the application of bias, however , np product departs from n, 2 and a net generation 
or recombination of carriers occur depending upcn the polarity. In case of a Schottky 
diode made on n-type silicon ,when the substrate is forward , biased electrons are 
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injected into the depletion region from the neutral bulk semiconductor and the holes 
are swept away from the metal . These additional electron-hole pairs combine in 
the depletion region to give a forward combination current . The recombination 
current can then be described by [4] 


4 = - 1] 

(1.15) 

qn.AdW 

with I ro = — 

2 r 

(1.16) 

and n t — \J . N c N v exp( £ ) 

(1.17) 


Here w is the width of the depletion region , r is the effective carrier lifetime 
within the depletion region , n t is the intrinsic carrier concentration , E g is the 
energy band gap and N c and N v are the effective conduction and valence band 
density of states, respectively. 

1.4.4 Minority Carrier Injection 

Under the low bias conditions , the Schottky barrier diode is a majority carrier 
device. But, at large forward bias , minority carrier current contribution increases 
and becomes more than the diffusion current. The electrons flow from the semicon- 
ductor into the metal under the forward bias and some of the holes are swept into 
the neutral region . The current I p caused by hole injection from metal into the 
neutral region of n-type semiconductor of donor concentration N<* is given by [ J/*] 


L = ^! lCTp( g)-i] 


NdL, 


(1.18) 


where D p is the, hole diffusion constant , L p = J. D p t p is the hole diffusion length 
and t p represents the hole lifetime in the neutral region. 
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Figure 1.8: Activation energy plots for determination of barrier height [7] 

1.5 Methods of Barrier Height Measurement 


In general , the barrier height is measured either using the current-voltage 
or the capacitance-voltage characteristics^ On the basis of thermionic emission- 
diffusion (TED) mechanism , the current through the Schottky barrier diode at 
a forward bias V is given by eqn (1.6). The plot of In (I) vs V is linear upto a 
point , for the forward bias V in excess of 3kT/q [4] . Afterwards the plot deviates 
showing saturation due to the series resistance term. In the plot of ln(I) vs V, 
the straight line portion is extrapolated to zero bias , so that the intercept at the 
ordinate is ln(I s ). By substituting the values of effective Richardson constant A** , 
saturation current I s in eqn (1. 5) , zero bias barrier height 4>b G can be determined 
at a temperature T for a given diode area Ad- The value of 4>bo is not very sesitive 
to the choice of A” since at room temperature a 100% increase in A** causes an 
increase of only 0.018Vin (j>b 0 [7]. 
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V (VOLTS) 

Figure 1.9: 1/C 2 versus applied voltage for determination of barrier height [7] 

Also , the activation energy (or the Arrhenius) plot i.e. In (I s /T 2 ) vs 1/T 
plot (Fig.1.8) can be made using saturation current data at various temperatures. 
(t>bo can be obtained from the slope of the best fit line [4]. Here the intercept at the 
ordinate is ln(A** A d ) where A d is the electrically active area. The main advantage 
in this approach is that the value of A d is not at all required for the evaluation of 
( fibo ■ 

Second method of determining the barrier height is by measurement of the 
capacitance of the Schottky diode as a function of reverse d.c. bias (V^) (Fig. 1.9). 
The charges of opposite signs are induced on the metal surface and in the depletion 
region of the semiconductor when a small a.c. voltage of few mV is superimposed 
on a d.c. bias. The relationship between the capacitance (C) and reverse bias Vr 
is given by[l] 
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C = A dt 


e s qN d 


\| 2 (fc - ^ -V R -f) 


(1.19) 


where e s is the permittivity of the semiconductor, N d is the dopant concentration 
, <j) n is the distance of the Fermi level from the bottom of the conduction band and 
other symbols have their usual meanings. Thus , if 1/C 2 is plotted against V R , a 
straight line with a slope of (2 / Aj 2 e s qNd) is obtained having an intercept on the 
voltage axis V R o which 1 /C 2 is zero. The value of barrier height can be than 
easily obtained as [4] 


4>b — Vro + 4 ( 1 - 20 ) 

q 

and <j) n , the depth of Fermi level below the conduction band , can be computed 
from the doping concentration N<*. can be evaluated from the slope of 1/C 2 vs 
Vi? plot too. 

1.6 Importance and Applications 

All the practical applications exploit the inherent high speed advantages of 
majority carrier conduction in Schottky diodes. The most important commercial 
application of Schottky diodes is in bipolar integerated circuits as clamps and to 
a lesser extent as load resistor substitutes , diode couplers and level shifters. They 
are also used as discrete microwave transistors , optical and nuclear particle de- 
tectors, etc. The Schottky diode is preferred over the pn junction diode because of 
high forward conductance per unit area , large reverse impedance , fast switch- 
ing speed about Ins , compact size ,and relatively easy fabrication technology at 
quite low temperatures. Metal-semiconductor rectifiers are relatively immune to 
the speed limitation associated with the minority carrier storage in p-n junction 
devices 7- ] . However, the reverse current in a Schottky barrier is quite large, 
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and is disadvantageous . Also , it is a surface device and hence sensitive to surface 
contamination and perimeter effects. 

1.7 Objective of the present work 

Extensive studies have been undertaken on different aspects of Schottky bar- 
rier diodes with the view to understand their characterstics better and find newer 
applications [3,4,'6'M). The most important aspect of metal-semiconductor con- 
tacts is the process which determines the flow of charges over the barrier from the 
semiconductor to the metal and vice- versa. Detailed knowledge of the conduction 
process involved is needed to extract the values of the barrier parameters. The I-V 
characteristics of Schottky diodes formed by depositing various metals such as Ir, 
Pt, Cr, Er, W ,Cu, and T i on silicon have been studied in detail [ 2.2-32.]. 

The fewer reports are available for Pd/Si Schottky diodes in particular and 
contain limited information on the conduction mechanisms ,e.g., Chin et. al. [$3] 
studied Pd 2 Si/Si Schottky barrier diodes on n-type silicon at 300K and on p-type 
silicon at 84 K . Werner and Guttler reported some data on Pd 2 Si Schottky 
barrier diodes while describing the temperature dependence of barrier heights for 
various metals on n-type silicon. Chand and Kumar have recently studied 

Pd 2 Si/n-Si contacts over a wide temperature range and interpreted the results on 
the basis of thermionic-emission-diffusion mechanism, assuming the presence of a 
Gaussian distribution of barrier heights at the interface. 

The current-voltage (I-V) characteristics of Schottky barriers measured at 
room temperature do not give adequate information about the conduction process 
and the nature of the barrier formed at the metal-semiconductor interface. Details 
of many possible effects that cause non-ideality in the diode I-V characteristics , 
and , in general , reduce the barrier height cannot be ascertained. The temperature 
dependence of the I-V measurements gives a better picture. Also, Schottky barriers 
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in devices operating in space and in some other applications are subjected to wide 
range of temperatures. So, information about their electrical characteristics over 
wide temperature range helps in understanding the actual A of these devices as well. 

An attempt has been made here to fabricate Pd/n-Si Schottky contacts 
and study their I-V characteristics in a wide temperature range (80K-300K). Since, 
fabrication involves deposition of palladium metal film and subsequent annealing , 
identification and growth of resulting phase(s) of silicide have also been explored 
using X-ray diffraction and reflectance measurements in UV-Visible range (200- 
900nm). In addition , simulation of (1) current-voltage (I-V) characteristics have 
been carried out for various mechanisms with a view to exmine their relative con- 
tributions and dominance , and, (2) apparent barrier height (BH) as a function of 
temperature . For the later case , barrier inhomogenities are believed to be present 
at the contact and described by a Gaussian distribution of barrier heights with a 
mean (fa ) and standard deviation (a 0 ). In this, uniqueness lies in considering the 
linear temperature dependence of the mean and variance (cr„). The purpose is to 
understand the effect of various parameters in detail. 



Chapter 2 
Experimental 

2.1 Fabrication of Schott ky diodes 

The Schottky diodes were fabricated upon n-type (phosphorus doped ) silicon 
wafers of (111) orientation . The n/n + silicon wafers had a 8-9 pm thick epitaxial 
n layer of resistivity ~2.6 ohm-cm(Nd = 3.9X10 15 cm -3 ) over the heavily doped 
n + region. The wafers were first degreased with soap solution and then were thor- 
oughly cleaned with organic solvents, viz, trichloroethylene , acetone and methanol 
in succession in an ultrasonic cleaner and rinsed in deionised water and dried. 

The back contact (ohmic) was made by depositing a layer of aluminium 
on the heavily doped n + side of the silicon wafer after usual cleaning and etching 
in dilute HF (HF:H2O=l:10) followed by washing in deionised water. A~ 200 nm 
thick aluminium deposition was carried out in a HHV coating unit model 12A 4D 

7 '/lew** 

under vacuum ~10~ 5 mbar, at 450° C for 30 minutes after introducing in a quartz 
tube , placed in a tubular furnace and connected to a pumping system model HHV 
VS114D. The silicon wafers were cleaned again from the epitaxial n or p layer 
(i.e. front) side in trichloroetylene , acetone and methanol and then etched in the 
dilute hydrofluoric acid (HF:H 2 O=l:10) and rinsed in deionised water for removal 
of silicon dioxide layer that usually forms on the silicon wafer during exposure to 
the ambient, aluminium deposition step and annealing. The back contact of the 
aluminium was protected from HF by laying a wax coating. Palladium film(~80 
nm thick) was deposited by employing an electron beam evaporation source in a 
HHV coating unit model 12A 4D-SC on to the cleaned epitaxial n-layer through 
circular holes of diameter 1 mm in a stainless steel mask. Subsequently, the samples 
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were taken out of the work chamber and put inside a quartz tube as before and 
annealed in vacuum of the order of~ 3-5 x 10 _6 mbar at 300^0, 400°C, 450°C, 500°C 
for various durations to form palladium silicide. The silicon wafers were finaly cut 
into small pieces such that each contained a palladium silicide dot and mounted 
on to a T05 header through silver epoxy. The contact to the top metal dot was 
established through a thin gold wire. 

2.2 Setup for I-V-T measurements 

An automated setup shown in figure 2.1 was used to record the I-V characterstics 
of Schottky diodes was over a wide temperature range (80K-300K). The bias was 
applied from a Keithley programable voltage source model 230 in steps of 5-lOmV 
and the current through the diode was recorded with Keithley picoameter model 485 
through a computer controlled system. To ensure steady state conditions, a delay 
of few seconds was observed between subsequent measuring steps. For recording 
the tempreture dependence of the I-V characterstics in the range 80K-300K a CTI- 
cryotronics close cycle helium refrigerator model 22 C equipped with a Lake Shore 
temprature controller model 330 was employed. 

2.3 Phase identification microstructure and re- 
flectance measurements 

For identification of the phase(s) formed and studying the growth kinetics 
of silicide, palladium film was deposited on a bare silicon wafer after usual clean- 
ing and etching in hydrofluoric acid exactly in manner described in section 2.1, 
but without mask. Samples were annealed at 200°C, 400°C, 450°C and 500°C for 
various durations in vacu um ^ 3-5xl0~ 6 mbar. Rich Siefert. ISO-Debyeflex 2002 X- 
ray diffractometer with CuKq radiation (A=1.5418JL) was employed to record the 
diffraction patterns. The calibration of the system was carried out with-a standard 
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Figure 2.1: Block diagr am of automated set-up for I-V measurements of Schott ky 
diode over a wide temperature range. 
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Picoammcter 
(Keithly 485) 

Figure 2.2: The circuit diagram for I-V measurements of a Schottky diode. 


silicon sample. The observed X-ray diffraction peaks were indexed by determining 
the interplanar spacings (d) corresponding to various planes using the well known 
Bragg’s law [<fd(2d sin 0 = X, 9 being the peak angle. The phase(s) identification 
was carried out by comparing the results with the standard data available in the 
powder diffraction filesj^i , 

The mi crostructure of the samples annealed at various temperatures for differ- 
ent lengths of time was examined in scanning electron microscope JEOL JSM 840A 
in the secondary electron(SE) mode. Also, reflectance measurments at near nor- 
mal incidence (by employing a 5° spacer) were carried out for the annealed samples 
in the wavelength range 200-900nm using a Hitachi UV 150-20 spectrophotometer 
equipped with a 60 <t> integrating sphere. 



Chapter 3 

Results and Discussions 

3.1 Forward I-V characteristics 

The forward I-V characteristics of Pd 2 Si Schottky diodes fabricated on n-type 
Si (111) as described in section 2.2 and annealed at 450°C for 45 min are shown in 
figure (3.1) at various temperatures in the range 80K-300K. Clearly, they are linear 
over several order of current. Further, they progressively become straight over a 
wide current range with decrease in temperature. Such an observation suggests the 
dominance of the thermionic emission -diffusion mechanism for current transport 
. A program in Basic is used to fit f he experimental I-V data in the thermionic 
emission-diffusion current equation (1.10). It involves iteration method and takes 
I s , 17, and R s as adjustable parameters. 

3.1.1 Zero bias barrier height , ideality factor and series 
resistance 

The zero bias barrier height (pbo) can now be evaluated from the I s data by 
substituting the values of the effective Richardson constant A**, diode area Aj, and 
temperature T. Thus , 4>bo is found at each temperature from equation (1.9). Figure 
(3.2) shows the zero bias barrier height as a function of temperature. It clearly 
depicts that q ^ decreases with fall in temperature steadily in the begining but 
quite rapidly below 120K.The variation of ideality factor and series resistance with 
the temperature are shown in figure (3.3) and (3.4) respectively . The ideality factor 
increases rather slowly with decrease in temperature but quite significantly at lower 
temperatures. There is a rapid increase of series resistance with fall in temperature 
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below 200K which is beleived to be due to factors responsible for increase in T) and 
/or lack of free charge carriers at lower temperatures [ 42 .]. 

3.1.2 Activation Energy Plot 

Equation (1.9) can be rewritten as 

^ (3.1) 

Therefore, the ln(I s /T 2 ) vs 1 /T ( i.e. activation energy or Arrhenius) plot should 
exhibit a straight line with slope giving the zero bias barrier height (cj>bo) and inter- 
cept at the ordinate determining the Richardson constant (A**) itself for the known 
diode area A d . This method gives a single value of barrier height. The experimental 
results , however correspond to the ln(I s /T 2 ) vs 1 /T plot , shown in figure (3.5). 
The data fits to a straight line only at high temperatures , i.e. , above 180K. The 
slope of the curve is infact varying continuously. 

3.2 Concept of Barrier Inhomogenities 

According to Chand and Kumar [ 36-3$ ] an abnormal decrease of barrier 
height (BH) and increase of the ideality factor ( 77 ) with the decrease in temperature 
cannot be explained adequately by incorporating interface states or interfacial oxide 
layer, tunneling , image force lowering and / or generation - recombination effects. 
A way to explain the abnormal behaviour of barrier height and ideality factor is 
to assume the existence of barrier inhomogenities and describe them with some 
suitable distribution function. Numerous evidences have since been reported for 
the existence of barrier height inhomogenities [ItZ-tty. Inhomogenities may be due 
to variation in thickness and composition of the silicide layer, nonuniformity of the 
interfacial charges and locally defective hot regions [-O]. 

Various types of distribution functions have been proposed to describe 
inhomogenities, e.g., a Gaussiah or a log iiormal. One can assume a Gaussian 
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Figure 3.1: Current - voltage characteristics of Pd/n-Si Sckottky contacts in the 
temperature range (a) 300K-180K (b) 165K-80K 
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Figure 3.2: Variation of the zero bias barrier height {4>bo) with temperature ,m the 
range 80K-300K 
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Figure 3.3: Temperature dependence of Ideality factor for Pd/n-Si Schottky con- 
tacts ,in the range &PK-300K. Notice the sharp increase in rj at low temperatures. 
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Figure 3.4: Temperature dependence of series resistance (R s ) for Pd/n-Si Schottky 
contacts ,m the range 80K-300K. Notice the sharp increase in R s at low tempera- 
tures. 



Figure 3.5: The ln(I s /T 2 ) vs 1000/T plots of Pd/n-Si Schottky contacts with best 
fit straight line ,in the range 80K-300K. It shows non-Unearity below 165K 
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Figure 3.6: The modified activation energy plots for a 0 0.1109 V and 0.063 V 
, together with the experimental data. The straight lines 1 and 2 indicate best 
fitting of data in temperature ranges 165-300K and 80-150K. 


distribution of barrier heights with a mean value of fa and standard deviation a , 
having the formfJ-3,*^,^ 

^ = ^ exp{JA ^ L) (a2) 

where preexponential term is the normalisation constant . The total current 
across the Schottky contact at a forward bias is then given by [43f 

I(V)= p I {faVMfaWb (3-3) 

J — oo 

where I (<j> b , V) is the. current for an elementary diode having barrier height <j> b 
at bias V and is given by equation (1.10) in case of thermionic emission-diffusion 
process. Substituting l((f>b , V) and p(<pb) from equations (1.10) and (3.2) in (3.3) 
and performing the integration from — oo to oo one gets 
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The barrier is known to depend on the electric field and hence on the applied 
bias [4]. Now , since barrier heights are represented by a Gaussian distribution, the 
entire profile should be affected by the applied bias. Hence considering the linear 
bias dependence of the distribution parameter, we have 

<}>b = 4>bo + 7 V and a = a 0 + £V (3.5) 

where 7 and £ are positive derivatives , and <j>^ and <j Q are mean barrier height 
and standard deviation at zero bias, respectively. Making these substitutions and 
neglecting the term involving £ 2 V 2 , equation (3.4) gets modified to 

1(V) = J.exp( 9( ^ fl » ) )[l - exp(- 9{V ~™') } (3.6) 

with 

I, = wr-expi-^.) 

Here I s is the observed saturation current , 0 ap and r) ap are the apparent barrier 
height and ideality factor respectively. <t> ap and r) ap are given by 

i>ap — 4>bo ~ 2 ~j7p an d ~ = (1 — 7 ) + —jjjT (3-7) 

It is clear that (j> ap depends upon the distribution parameters (i.e. mean 6b 0 and 
standard deviation o 0 ) and temperature T. Obviously, decrease in the barrier height 
is caused by mere existence of distribution and its effect becomes more and more 
pronounced at low temperatures. On the other hand, is a function of rate of 
change of both the mean 4>bo and the standard deviation <r 0 with the bias besides cr 0 
itself. 


3.3 Interpretation of the Experimental Data 


The diode parameters derived from the fitting of the experimental data in 

$ 7 - \ h 

equation (3.6) can therefore be interpreted using equation (3.1) and (3.£) . Thus 
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the <j> ap vs 1/T plot should be a straight line giving ^ and a 0 from the intercept 
at the ordinate and the slope respectively. The I-V data however, reveal , two 
straight lines in <j> ap vs 1/T plot giving two values of (0.93V and 0.63V) and a Q 
(0.118 and 0.067V). Such an observation is , in fact attributed to the presence of 
two Gaussian distributions of barrier height at the interface and discussed in length 
later in section 3.4 

Also combining equations (3.6) and (3.1) one obtains 

Mj |) - $*5 = WA.A") - ^ (3.8) 

The nonlinearity in the activation energy plot can now be understood on the 
basis of equation (3. 8 )■ The factor ln(|t) — is calculated to modify the experi- 
mental data and plotted as a function of 1/T to obtain a straight line corresponding 
to a single activation energy around accurately. The intercept of the straight 
line at the ordinate gives the Richardson constant itself. Figure (3.6) shows the 
modified activation energy plot . Notice that the deviation is removed , but, the 
data fits with two straight lines. These correspond to the value of ^ 0.93 and 
0.63V and Richardson constants 0.99 x 10 6 and 1.69 x 10 6 Am ' 2 K~ 2 . Similarly, the 
temperature dependence of ideality factor (r} ap ) can be understood on the basis of 
equation (3.7). It indicates that l/rj ap vs 1/T plot should yeild a straight line with 
slope giving <r 0 q£/k and intercept (1-q). The increase of rj with fall in temperature 
arises due to the term containing £ (i.e., the rate of change of a with bias) . The 
plot of l/r) ap vs 1/T also contains two straight lines . Both correspond to negative 
slopes making £ essentially negative; the values being -2.38mV and-0.14mV in the 
temperature range 80-135K and 150-300K, respectively. 

It may be mentioned that the studies on many other simultaneously fabri- 
cated diodes revealed slight variation in the values of the mean barrier height and 
standard deviation at zero bias. Such observations are primarily attributed to the 
local chemistry and conditions of the interface. 
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3.4 Evidence for a Double Distribution 

Figure (3.7) shows yet another lnl vs V characteristics of a Schottky diode fabri- 
cated after annealing at 400°C for lh in vacuum. These clearly depict increasingly 
linear behaviour over several orders of current and shift towards higher bias side 
with decrease in temperature . The variation of apparent barrier height and ideality 
factor with temperature as obtained in this case are given in figures (3.8) and (3.9) 
respectively. According to equation (3.J) , the plot of 4> ap should be a straight fine 
with intercept at the ordinate determining the zero bias mean barrier height <j>bo 
and the slope giving the zero bias standard deviation a 0 . However, when the <p ap 
vs 1/T plot is made (figure (3.10)) , data fits nicely , as in the previous case , with 
two straight fines (transition occurring ~ 135K), implying thereby that <f> bo and a 0 
assume different values in the temperature range 165-300K and below 165K. The 
intercepts and slopes of these straight lines yield values of 4 ^ and a 0 as 0.93V and 
0.112V in 165-300K and 0.629V and 0.067V below 165K. 4> ap can now be estimated 
from equation (3.7) itself over entire temperature range 80-300K for each set of 4>bo 
and cr 0 values. This is presented in figure (3.11) by continuous curves. 

Similarly , 1 /r) ap vs 1/T plot should also possess different characterstics 
in the two temperature ranges if contact contains two barrier heights distribution. 
This is indeed the case as data fits with two straight lines in the temperature ranges 
300-165K and below 165K ( see ,e.g. , figure (3.12)). Their slopes yield qcr 0 £/k and 
intercepts at the ordinates give (1-7 ) . The value of £ extracted from the slopes 
are -2.4mV and -0.14mV whereas the values of 7 obtained from the intercepts are 
-0.75 and 0.187. 

Figure (3.13) shows the conventional activation energy ln(I s /T 2 ) vs 1/T 
plot . The experimental data are seen to fit asymptotically with a straight fine 
at higher temperatures only. The zero bias barrier height obtained from the slope 
of this straight fine comes out to be 0.34V . Likewise the value of Richardson 
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Figure 3.8: Variation of the zero bias barrier height with temperature ,in the 
range 80-300K. 
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Figure 3.9: Temperature dependence of Ideality factor for Pd/n-Si Schottky con- 
tacts ,in the range 80-300Jt. Notice the sharp increase in rj at low temperatures. 
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Figure 3.10: The barrier height (of, 0 ) obtained from I-V measurements as a function 
of inverse temperature.Notice the transition occuring around 150K . 



Figure 3.11: The zero-bias barrier height as a function of temperature.Full curves 1 
and 2 correspond to <j> ap calculated using equation (3.8)for two distributions having 
fa (0.93V and 0.629V ) and <r 0 (0.112V and 0.0667V}. 
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Figure 3.12: The inverse ideality factor (1/^) versus inverse temperature (1/T) 
plot. The data show linear variation in the two temperature ranges with transition 
around 200K 



Figure 3.13: Conventional activation energy ln(I s /T 2 ) vs 1000/T plot. Notice the 
deviation from linearity below 180K. 




1000/T (K -1 ) 

Figure 3.14: The modified activation energy plots for <r D 0.112 V and 0.067 V 
, together with the experimental data. The straight lines 1 and 2 indicate best 
fitting of data in temperature ranges 165-300K and 80-150K , respectively. 

constant determined from the intercept at the ordinate is 3.34x10 Am~ 2 K~ 2 . Both 
these parameters are quite unreasonably low as the reported value of the barrier 
height for Pd 2 Si/nSi Schottky diodes is 0.735V [<$£] and Richardson constant for n- 
type silicon is 1.12xl0 6 Am ~ 2 K~ 2 [4] . Moreover, the activation energy plot shows 
nonlinearity at low temperatures. To explain these discrepancies the I-V data are 
interpreted in terms of equation (3.10) and a modified activation energy plot (figure 
(3.14)) is made. For this the term ln( j5?) — is calculated for both the values of 
cr 0 ( associated with the Gaussian distributions of barrier heights) found above and 
then plotted as a function of inverse temperature . The best linear fitting yields 
zero bias mean barrier heights as 0.969V and 0.680V which match exactly with 
the mean barrier heights obtained from the 6 ap vs 1/T plots earlier. The intercept 
at the ordinate is just ln(A<j A**) and gives Richardson constant A** as 1.22xl0 6 
Am ~ 2 K~ 2 and 1.69xl0 6 Am -2 K~ 2 for the two straight lines drawn. 

All the considerations above suggest that the obsersved experimental data can 
be satifactorily explained by assuming the existence of two Gaussian distributions 
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of barrier heights in the present Pd 2 Si/nSi Schottky diodes. 

3.5 X-ray Diffraction Studies 

Fig. (3. 15), (3.16), (3. 17), (3. 18), and (3.19) show X-ray diffractograms of grow- 
ing silicide at different stages of annealing at 200°C , 400°C ,450 o C and 500°C . 
These are recorded with CuKa radiation(wavelength A = 1.5418 A ) at the same 
sensitivity and time constant. The Bragg angles (26) , d- values and relative inten- 
sities of various reflections are summarized in Table (3.1) and (3.2) for palladium 
and palladium silicide , respectively. Palladium film corresponds to f.c.c. phase 
with a=3.89A ^}]. The diffraction data suggests that annealing of Pd/n-Si<lll> 
at 400-500*0 invariably leads to silicide of composition Pd 2 Si having Fe 2 P struc- 
ture with hexagonal unit cell ; the parameters being a=13.055A , c=27.490A $t}. 
Also , there is a tendency of preferred orientation , e.g. , at 500°C , Pd 2 Si grains 
assume preferred orientation such that their basal plane lies parallel to (111) of 
silicon substrate. (Notice appearance of 00.16 reflection and progressive increase in 
its intensity with annealing) . This means <111> Si is parallel to <00.1> Pd 2 Si. 
This result is consistent with the findings of Buckley and Moss [i'-J. Further, at 
400°C and 450°C , the intensities do not match with the ASTM data |%t| of poly- 
crystalline Pd 2 Si and therefore indicate the presence of preferred orientation. The 
diffraction data suggests that (22.8) Pd 2 Si lie parallel to (111) Si. With lapse of 
time , some readjustment occurs and grains assume other orientations as well and 
move towards true polycrystallinity ,i.e., Pd 2 Si exhibits epitaxy in the begining but 
turns into polycrystalline disposition as it thickens. 

Palladium reacts with silicon in accordance with the equation 2Pd + Si — > 
Pd 2 Si. Thickness d* of Pd 2 Si formed can be estimated from 

j* Mp d ,siPpd 

dpd ~ AnNp&pw * 1 





4i 



Figure 3.16: Diffractogram for the sample annealed at 200°C for differ- 
ent lengths of time. 
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Table 3.1: Diffraction peaks for Pd 


Pd 

Reflections 

hkl 






111 

200 

220 

311 

222 

as deposited 

29 

40.1 

47.1 

68.3 

82.5 

86.6 


dO) 

2.248 

1.930 

1.37 

1.171 

1.122 


R(I) 

85 

15 

30 

10 

8 

annealed 

29 

40.1 

47.1 

68.3 

82.5 

86.6 

at 200° C 


2.248 

1.930 

1.37 

1.171 

1.122 

lh 

R(I) 

80 

20 

40 

15 

5 


(2 9 

40.1 

47.1 

68.3 

82.5 

86.4 


dO) 

2.248 

1.930 

1.37 

1.171 

1.126 

lOh 

R(I) 

80 

18 

35 

10 

3 

annealed 

29 

40.1 

47 

68.2 

82.4 

86.4 

at 400°C 

dO) 

2.248 

1.933 

1.375 

1.173 

1.126 

0.5h 

R(I) 

180 

30 

60 

20 

10 


29 

40.1 

47 

- 

82.4 

_ 


dC4) 

2.248 

1.933 

- 

1.173 

— 

6h 

R(I) 

80 

15 

- 

10 

- 


29 

40.1 

- 

- 

82.3 

_ 


d{A) 

2.248 

- 

- 

1.172 

— 

lOh 

R(I) 

30 

- 

- 

5 

- 


29 

40.1 

- 

- 

- 

_ 


dC4) 

2.248 

- 

- 

_ 

_ 

16h 

R(I) 

20 


- 

- 

- 

annealed 

29 

40.1 

- 

68.2 

82.5 

- 

at 450°C 

<&) 

2.248 

- 

1.375 

1.171 

_ 

lh 

R(i) 

50 

- 

28 

8 

- 


29 

40.2 

- 

68.3 

82.1 

- 


d(&) 

2.246 

- 

1.373 

1.174 

- 

3h 

R(I) 

55 

- 

26 

7 

- 


29 

40.1 

- 

68.1 

82.5 

- 


diA) 

2.248 

- 

1.377 

1.171 

_ 

5h 

R(I) 

40 

- 

20 

5 

- 


29 

40.2 

- 

- 

- 

- 


dfl) 

2.246 

- 

- 

- 

_ 

7h 

R(I) 

20 

- 

- 

- 

- 

annealed 

20 

40.1 


68.5 

82.3 

- 

at 500°C 


2.248 

- 

1.369 

1.172 

- 

0.5h 

R(I) 

50 

- 

30 

10 

- 


29 

40.1 

- 

68.3 

82.5 

- 


d<& 

2.248 

- 

1.373 

1.171 


1.5h 

R(i) 

44 

- 

28 

7 

- 


29 

40.2 

- 

68.8 

82.4 

- 


<mA) 

2.246 

- 

1.364 

1.1715 

- 

2.5h 

R(I) 

39 

- 

24 

5 

- 


29 

40.1 

- 

68.2 

- 

- 


cKA) 

2.248 

- 

1.375 

- 

- 

5h 

R(I) 

26 

- 

20 

- 
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Table 3.2: Diffraction peaks for Pd 2 Si 


Pd 2 Si 

Reflections , hkl 



228 

408 

0016 

221j§ 

2416 

6016 

as deposited 

26 

- 

- 

- 

- 

- 

- 



- 

- 

- 

- 

- 

- 


R(i) 

- 

- 

- 

* 

- 

- 

annealed 

26 

- 

- 

- 

- 

- 

- 

at 200° C 

d$) 

- 

- 

- 

- 

- 

- 

lOh 

R(I) 

- 

- 

- 

- 

- 

- 

annealed 

26 

38.1 

- 

- 

- 

- 

- 

at 400° C 

d(A) 

2.362 

- 

- 

- 

- 

- 

0.5h 

R(I) 

30 


- 

- 

- 

- 


26 

38.1 

- 

53.7 

- 

- 

74.8 


dO) 

2.362 

- 

1.709 

- 

- 

1.270 

6h 

R(I) 

140 

- 

9 

- 

- 

24 


26 

38 

41.5 

53.8 

- 

70.3 

74.8 


d(A) 

2.366 

2.176 

1.703 

- 

1.339 

1.270 

lOh 

R(I) 

180 

32 

10 

- 

40 

26 


26 

38.1 

41.4 

55.5 

61 

70.4 

74.9 


d(l) 

2.362 

2.181 

1.713 

1.521 

1.337 

1.268 

16h 

R(I) 

200 

40 

12 

8 

50 

30 

annealed 

26 

38 

- 

53.3 

- 

70.3 

- 

at 450°C 

m 

2.366 

- 

1.718 

- 

1.339 

- 

lh 

R(I) 

25 

- 

10 

- 

24 



26 

38.1 

- 

53.3 

- 

70.3 

74.9 


<m 

2.362 

- 

1.718 

- 

1.339 

1.268 

3h 

R(I) 

35 

- 

12 

- 

23 

5 


26 

38 


53.4 

- 

70.4 

74.8 


d&) 

2.366 

- 

1.715 

- 

1.337 

1.270 

5h 

R(I) 

40 

- 

11 

- 

23 

10 


26 

38.1 

- 

53.3 

- 

70.3 

74.9 


d&) 

2.362 

- 

1.718 

- 

1.339 

1.268 

7h 

R(I) 

60 

- 

12 

- 

30 

16 

annealed 

26 

- 

- 

53.7 

- 

70.3 

- 

at 500° C 

d&> 

- 


1.706 


1.339 

- 

0.5h 

R(I) 

• 

- 

25 

- 

18 

- 


26 

38 

- 

53.4 

- 

70.4 

- 


d&) 

2.366 

- 

1.715 

- 

1.337 

- 

1.5h 

R(I) 

20 


38 

- 

19 

- 


26 

38.1 


53.5 

60.9 

70.4 

- 



2.362 

- 

1.713 

1.521 

1.337 

- 

2.5h 

R(I) 

24 

- 

42 

14 

19 

- 


26 

38.1 

41.5 

53.3 

60.9 

70.3 

74.9 


d&> 

2.362 

2.176 

1.718 

1.521 

1.339 

1.268 

5h 

R(I) 

37 

5 

47 

15 

18 

10 
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where Mp^si and A p d are the molecular weight of Pd 2 Si and atomic weight of 
Pd , respectively, p Pd and pp d2 s% their densities , d“ Pd thickness of Pd-film is deter- 
mined by optical interferometry as ~1050 ± 50A. Substituting the values of various 
paramaters as M Pd2Sx =240.89 , A Pd =106.4 , N Pd = 2 , p Pd = 12.02g/cm 3 p Pd2 s t = 
9.462g/cm 3 , d p d2S% is estimated as 1510 ± 70 A. 

Wittmer and Tu \0\ have shown that the growth of Pd 2 Si on silicon wafer 
of <100> and < 11 1 > orientations is diffusion limited with activation energy (E a ) 
lying in the range 1.3-1. 5 eV. Also, Pd 2 Si follows parabolic growth law(i.e., d** = 
2(Dt) 1</2 , where d~ is the silicide thickness, t is the annealing time and D(T)=D 0 
exp (-E a /kT) is the chemical interdiffusion constant; the value of D 0 is 0.48 or 
0.18 cm 2 /sec for p-doped Si <100> or <111> orientation) and is independent of 
the substrate orientation.The rate limiting step is interdiffusion through the silicide 
layer rather than the interfacial reaction . Further, the dominant diffusing species 
is palladium in Pd 2 Si via vacancies. 

These observations are at variance with the earlier reports of Cheung and 
Nicolet [51] who found diffusion constants to be in the range 10~ 15 - 10 -13 cm 2 /s 
at 240 o C and activation energy of Pd 2 Si formation as 0.9i0.1eV. Our value of 
interfacial diffusion constant lie in the range 2. 2-6. 3 x 1CT 15 cm 2 /s at 400-500 c 'C and 
are closer to the later case. 

3.6 Reflectance measurements and SEM studies 

Fig. (3.20) , (3.21), (3.22) show reflectance(%R ) of Pd/n-Si system as such 
and after annealing for different lengths of time at 400, 450 and 500°C, respectively. 
Accordingly at 400°C, the reflectance spectra preserve the trend but decreases with 
time of annealing . On careful examination, it was found that .patches of Pd 2 Si 
developed first at a number of places , they enlarged to islands and eventually 
covered the entire area by coalescence. At 450 °C , the reflectance (%R) decreases 


Reflectance(%) 
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Wavelength(nm) 


Figure 3.20: Reflectance versus wavelength plots for the sample annealed at 400°C 
for different lengths of time 




Wavelength(nm) 


Figure 3.21: Reflectance versus wavelength plots for the sample annealed at 450°C 
for different lengths of time 


Reflectance(^) 



Wavelength(nm) 


Figure 3.22: Reflectance versus wavelength plots for the sample annealed at 500°C 
for different lengths of time 
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initially , then increases upto 2h and eventually decreases at all wavelengths. In 
this case no patches , as observed at 400 °C , were found . Instead it is believed that 
lateral growth of Pd 2 Si occurs. While the first decrease can be attributed to initial 
growth of Pd 2 Si , the increase in the reflectance may be due to recrystallization and 
emergence of preferred orientations. As the time lapses , Pd 2 Si exhibits increased 
polycrystallinity (section 3.5) and so reflectance decreases as usual. At 500 °C , 
reflectance decreases with annealing time as basal plane (00.1) is predominantly 
lying parallel to the underlying (111) plane of silicon . But , after 1.5h , there 
is a slight increase in reflectance presumably due to other orientation developing. 
Eventually randomization occurs and we observe decrease of reflectance after 2.5h. 

The scanning electron micrographs are shown in fig(3.23) 


CENTRA l'SRAHI 

1 1 T XASSIpy# 


4mt+ A 


1 





Results and Discussions 


52 





Figure 3 23' Scanning electron micrograpbsof Pd/n-Si Schottky contacts showir« 
STstStnre (a)after annealing at 400»C for 12h (b) after annealmg at «0»C for 
4h (c) after annealing at 500*C for 30 min. 


Chapter 4 
Simulation Studies 

4,1 I-V characteristics 

In order to understand the relative contributions of various current mechanisms, 
in Sehottky diodes, I-V characteristics have been simulated for all of them using 
known parameters. For this, thermionic emission-diffusion, generation-recombination 
and tunneling processes are chosen. The thermionic emission-diffusion current through 
a Sehottky diode under a forward bias of V volts is given by [4,7]. 

he = IJgXp (|^) ~ 1] (4.1) 

with saturation current I s given by 

I, = AtA-T^pf^) (4.2) 

where A** is the Richardson constant (1.12 x lOP Am~ 2 K~ 2 ) and o b is the barrier 
height (for PdzSi/n — Si contact, its value is ~ 0.735eV [11]) . is diode area 
(for 1mm size dot Ad = 7.85x10 _7 tu 2 ) and k the Boltzmann constant is 1.38 x 
10 -23 JjK . With these parameters, I s has been simulated using eqn(4.2) at various 
temperatures in the range 50K-300K. For each I s , in turn , thermionic emission- 
diffusion current (I te ) can be computed from eqn(l.lO) at various forward bias levels 
at the corresponding temperature. Figure 4.1 shows the results ln(I) vs (V) plots 
at various temperatures. Notice that saturation current I s lies in the region ~ 
10 -36 .A — 10~ 9 A at temperatures 100K - 300K. 

The generation-recombination current through a Sehottky diode under for- 
ward bias of V volt is given by [4]. 
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Figure 4.1: Simulated la(I) vs V plots for the current due to thermionic emis- 
sion-diffusion process at different temperatures 
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, / 

Ir = -Mezp^} - 1] (4-3) 

with saturation current current 

1,0 = qn l wA d /2r (4.4) 

where, n, is the intrinsic carrier concentration, w is the width of the depletion region, 
A d is the area of the junction(interface) andr is the carrier life time, n, is given by 

[7] 


n. = (4.5) 

wh('re in* and m£ are the effective mass of the electron and hole respectively 
that have value 1.18m and 0.81m ; la being the mass of electron , E g is the band 
gap of silicon, k the Boltzmann constant and h is the Plank’s constant . The value 
of n, is determined at different temperatures with E g = 1.12eV , w = 0.3/im, Ad — 
7.85xl0" 7 m 2 and r = 10 _6 s [7]. Generation recombination current can be obtained 
from eqn(4.3) and eqn(4.4) at different temperatures (Fig.4.2). Notice that the 
saturation current is always 2-5 orders of magnitude less than the corresponding 

values obtained for the TED process (Fig. 4.1). 

The tunneling current through the barrier under a forward bias of V volt is 

given by [4, |X,f3] 


It = It 0 exp(V/E') 


(4.6) 


with the saturation current Ito 

Ito ~ 

and every parameter 




k ^ c° sh21 ffl 


(4.7) 




(4.8) 
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Figure 4.2: S imul ated In (I) vs V plots for current due to generation recombination 
process at different temperatures 
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Eo and arc related such that 


E 0 = Eoocoth 


r lE M 

~kT 


(4.9) 


withE^ given by 


An m*€ s 


(4.10) 


where N r ; is the donar density and e s is permittivity of silicon. At donar con- 
centration Nd ~ 1.35 x 10 21 ,m* = 1.18m and e s = 11.9, E C<J comes out to be 
1.88 x 10 4 eV. Using this value of Eq^E' can be calculated at different temper- 
atures and in turn,/ to with A * — 1.12 x 10 6 AK~ 2 m 2 ,6f, = 0.763eV\ Thus, tunneling 


current is obtained at different bias from eqn(4.6) for temperature range 50K-300K . 
Fig(4.3) shows variation of tunneling curent with bias at various temperatures. The 
plots reveal very little variation of current with bias. Also, the current level drops 
down significantly with decrease in temperature. Fig (4.4) shows bias dependance 
of current more clearly at 300K. Obviously, tunneling current is always very low 
in comparison to current resulting due to TED and generation recombination pro- 
cesses.(See eg., Fig (4.5)). 

The ln(I)-V plots should be linear for all the three currents as per the govern- 
ing equations (4.1), (4.3), (4. 6). It has however been observed that the experimental 

o • 1 /, «- r 'c 

ln(I)-V plots are linear over several orders of current^higher biases.Such an effect is 
actually caused by the series resistance ( R s ) present .In this situation , bias applied 
to the diode gets modified to (V-LR s )in equations(4.1)(4.3)(4.6) Saturation effects 
appear and the current can be estimated by iteration only as described in chapter 
3. Also, fitting of I-V data in TED current equation gives barrier height 4> ap and 77.lt 
is usually found that <f> ao decreases and r) increases with decrease in temperature 
and activation energy plot deviates from linearity at low temperatures.lt amounts to 
saying that the actual current is more than the contribution made by TED process 
itself. These anomolies are at first sight indicate a deviation from the thermionic 
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Figure 4.3: Simulated ln(I) vs V plots for the current due to tunneling process at 
different temperatures 
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emission-diffusion mechanism.The possible reason could be the simultaneous oper- 
ation of other processes viz., tunneling through the barrier or recombination in the 
depletion region. 

The tunneling involves field emission and thermionic field emission through 
the barrier. The process applicable, however, depends upon the parameter 


Voo = ^-(^) 1/2 = 18.5 x 10 15 (^-)V 2 eF (4.H) 

4tt m* e m e e r v ' 

The field emission becomes important when E !Xj kT /q, whereas thermionic 
field emission (TFE) dominates when E^ ~ kT / q and thermionic emission-diffusion 
if Eoo < kT/q [4], In the present case E^ = 0.188mV with N D = 1.35xl0 21 m~ 3 
, m s = 1.1 and e s = 11.9.This leads to qEoo/kT equal to 0.0071 at 300K and 0.027 at 
80K . Obviously , the condition prevailing in both the situations is E m ;§> kT/q(~ 
0.259 mV and ~ 0.7 mV at 300K and 80K, respectively) and hence the possibility of 
FE and TFE can be easily ruled out. The tunneling current as estimated from eqn. 
comes out to be very low ~ 10 _21 A at 300K and ~ 10 _47 A at 100K, whereas, the 
experimentally obtained values just at forward bias of 0.1V are 10 _6 A and 10 -11 A 
at 300K and 100K respectively. 

The second mechanism , i.e. , the generation and recombination requires (see 
,e.g., eq(4.3)) the ideality factor to be close to 2. But this value of q is found at 80K. 
Below 80K,the ?? remains much smaller than 2.The saturation current calculated 
using eqn(4.4),comes out to be ~ 10" 12 Amp at 300K and 10~ 30 A at 100K. Obviously 
the estimated current values for recombination process are excessively low compared 
to those actually observed in experiments. Chen et al [52] have pointed out that 
if recombination current exists, the barrier height evaluation on the basis of TED 
theory yields progressively lower values with decrease in temperature;for example, if 
I r constitute 91% of the total current at (say) 100K, barrier decrease expected is 
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* v = &(0)-|2+a 1 T-2^MI (4.15) 

Thus, (j) aj) contain terms amounting to direct and inverse temperature dependence. 
The a oove equation is used to determine o ap in the temperature range 5 OK- 35 OK 
for different values of ct x and « 2 . For this data taken are 4,(0 K) = 0.8V and <r o (0) 
= 0.06V. The value of a x is varied from -0.3mV/K to +0.3mV/K and that of a 2 
from -3(mV) 2 /K to +3(mV) 2 /K. 

Notice the coefficient a 2 is for variance cr 2 . The corresponding values 
of temperature coefficient of standard deviation lie in the range -0.025mV/K to 
+0.25mV/K . Fig. (4.5) shows the variation of apparent barrier height as a function 
of temperature at various values c*} when a 2 is equal to zero. The situation with 
«i = 0 corresponds to equation(4.12) itself. It depicts decrease of barrier height 
with decrease in temperature due to the presence of Gaussian distribution of barrier 
heights with mean(^»j, 0 ) and variance (cr 2 ) . When itself varies with temperature 
(i.e., CV] = -3 to 3 mV/K and a 2 = 0 ) the extent of decrease is reduced or enhanced 
depending upon the sign of temperature coefficient (i.e., positive or negative). The 
difference is , however, more prominent at higher temperatures ,i.e., above 150K. 

The effect of a 2 on the barrier height is just opposite and shown in Figs. (4.6) 
and (4.7) for Qi = -3mV/K and 3mV/K, respectively. Here , the extent of decrease 
is less or more if sign of a 2 is negative or positive. Physically, it means that when 
a 2 is negative (equivalent to narrowing of the distribution function, i.e., moving 
towards homogenization of barrier heights ) , the decrease of barrier height will 
be less than predicted by equation(4.12). On the other hand , when a 2 is positive 
distribution function widens and amounts to increase in the value of <r 0 with 
rise in temperature . Consequently, there is effectively more reduction in <f>bo then 
determined by equation(4.12). 
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Figure 4.6: Apparent barrier height versus temperature plot with varying from 
-3 to 3 mV /K and a 2 zero 
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Figure 4.7: Apparent barrier height versus temperature plot with a#-3 mV/K and 
ct 2 varying from -3 to 3 mV 2 /K 
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Figure 4.8: Apparent barrier height versus temperature plot with a a :3 mV/K and 
a 2 varying from -3 to 3 mV 2 /K 



Chapter 5 
Conclusions 


• Pd/n-Si(lll) Schottky barrier diodes can be prepared successfully by depositing 
a palladium layer of ~ 80CL4 thickness using an electron beam evaporation source 
in vacuum ~ 10 -6 mbar onto a thorougly cleaned silicon substrate. 

• I-V characteristics of Pd/n-Si(lll) Schottky contacts can be explained on the 
basis of thermionic emission-diffusion (TED) mechanism . The ln(I) vs V plots 
over the temperature range (80K-300K) exhibit linearity over several orders of cur- 
rent and saturates due to effectiveness of voltage drop occuring across the series 
resistance (R s ) of the diode. 

• Several diode parameters , viz, barrier height , ideality factor , series resistance 
and saturation current can be determined by the fitting of the I-V data in the TED 
current equation using a computer program. While the zero bias barrier height (by, 
decreases , ideality factor rj and series resistance R s increase with the decrease in 
temperature. Moreover , the changes in <f> bo , rj , and R s become quite significant at 
low temperatures. 

• The abnormal behaviour of Pd/n-Si Schottky diodes i.e. , decrease of zero 
bias barrier height and increase of ideality factor with decrease in temperature can 
be successfully explained on the basis of TED mechanism by incorporating the 
concept of barrier inhomogenities and describing them by a Gaussian distribution 
of barrier heights. The values of parameters obtained are : zero bias mean barrier 
height (<j)bo) 0.93V and 0.63V , standard deviation ( a Q ) 0.112 and 0.067V in the 
temperature range (165K-300K) and (80K-165K), respectively. 

• The <t>b 0 vs 1/T plot yields straight line(s) to demonstrate the number of 
barrier distribution(s) , viz, single, double or multiple existing at the Schottky 
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contact giving values of both and o 0 . Similarly, inverse ideality factor (1/ 77) vs 
inverse temperature (1/T) plot displays linearity and it’s intercept at the ordinate 
and slope determine the bias coefficients (7 and £) of the mean barrier height and 
the standard deviation , respectively. 

• The annealing of Pd/n-Si <11 1> at 400-500°C invariably leads to silicide 
formation of the composition Pd 2 Si having a hexagonal structure (Fe 2 P type) with 
parameters a = 13.055A , c= 27.490JL. Also, there is a tendency of preferred 
orientation in the begining but with lapse of time , some readjustment occurs and 
grains assume other orientations as well to exhibit true polycrystallinity. 

• Simulation studies suggest dominance of TED mechanism in the Pd/n-Si 
Schottky contacts. Also, negative temperature coefficients of mean barrier height 

and variance (<r^) of the Gaussian distribution function (representing the bar- 
rier inhomogeneties at the Schottky contact ) are found to have opposite effects on 
the temperature dependence of apparent barrier height. While the former causes 
decrease in barrier height the later compensates the change by a fixed amount. 
The overall effect is that Gaussian distribution function shifts towards lower bar- 
rier height and at the same time becomes narrower with rise in temperature , i.e., 
effective reduction in barrier height occurs. 
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